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1.0 INTRODUCTION 


Obtaining global knowledge of the earth’s gravitational fine structure is 
necessary for understanding solid earth and ocean dynamics, and is a major 
objective of NASA's Earth and Ocean Dynamics Application Program [13]. Several 
satellite mission configurations and data types have been proposed for this task, 
including satellite altimetry, satellite gravity gradiometry, and satellite-to- 
satellite tracking. A review of these concepts for gravity field determination 
has been given by Argentiero and Lowrey [3]. The feasibility of utilizing a 
satellite-to-satellite tracking mission configured as two low altitude drag-free 
GRAVSAT spacecraft in identical orbits ("Low-Low") has been investigated by 
Estes and Lancaster [7]. This publication will be refered to as report [A]. The 
present study will employ the error analysis techniques of report [A] to investigate 
the feasibility of utilizing a "high-low" SST mission involving a high altitude 
relay satellite and a low altitude GRAVSAT spacecraft to recover the global 
gravitational fine structure. 

The 10 cm accuracy requirement placed on the gravity fine structure (see [13]) 
necessitates an extremely large number of parameters in the field model. A 
simultaneous data reduction for all model parameters would be a tremendously 
difficult computer task, even with the most sophisticated new computer systems. 

This supports the idea of representing the gravitational fine structure with a 
local mathematical model and estimating a local area of gravitational perturbations 
using local data. This would, of course, require a fewer number of model parameters 
and make the problem more manageable. The global solution could then be built 
up from a set of these local solutions. This idea has been tested by simulation 
for SST data in report [A] and by Schwarz [15] and Hajela [8]. However, while 
adjusting a set of the model parameters for a local recovery, the remaining 
unadjusted parameters must be constrained to a priori values. Unless the gravity 
field parameterization and the data exhibit "orthogonality", the uncertainties, 
of the unadjusted parameters will corrupt, or "alias", the estimates of the 
adjusted parameters. Orthogonality has been defined by Argentiero [2] as follows: 
For a given estimation procedure, the j th adjusted parameter is orthogonal to 
the kth unadjusted parameter if the aliasing contribution of k to j is zero. 

The mathematics of this concept is described in Section 4.0. 
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The concept of the "high-low" GRAVSAT SST mission is that the low altitude 
satellite is strongly perturbed by the gravitational fine structure while the 
high satellite is virtually unaffected, resulting in the data obtained from 
tracking the relative motion between them being highly sensitive to geopotential 
variations. Only a few geopotential terms are important at the 10 cm level for 
altitudes above about 5 earth radii (for non-resonance conditions). 

The low altitude (250-350 km) GRAVSAT is a drag-free spacecraft constrained 
to follow a purely gravitational path in a polar circular orbit. The orbit 
geometry provides global data coverage. The high altitude satellites considered 
in this study are the GEOPAUSE [14] and geosynchronous. 

The purpose of this investigation is to apply covariance error analysis 
techniques to determine under what estimation strategies the "high-low" mission 
is feasible for determining the short wavelength gravity field uniquely and 
accurately from local recovery of gravitational fine structure when the effects 
of aliasing and satellite state uncertainties are considered. A multi-short 
data arc approach to the problem is utilized (see report [A], and [4]) whereby 
the satellite epoch states are estimated together with the gravity parameters. 

For local solutions, long arcs lead to severe aliasing problems. The results 
obtained for the "high-low" configuration are compared with "low-low" results 
from report [A] in Section 6.2. To help establish the feasibility for obtaining 
global solutions from a composite of local solutions, a set of simulated global, 
long arc solutions for the "low-low" configuration are presented in Appendix A, 
and show compatibility with the short arc, local "low-low" results of report [A]. 
The drag compensation system of the GRAVSAT satellite has been assumed error 
free in this study, so that the results are somewhat optimistic. In a more 
detailed study the aliasing effect in the parameter estimation due to imperfect 
knowledge of the surface force compensation must be included. 
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2.0 SURFACE DENSITY GEOPOTENTIAL REPRESENTATION 


The purpose for the present study of the high-low SST experiment is to determine 
the feasibility of recovering the global fine structure of the earth's gravity 
field as a composite of recoveries over local regions. The total geopotential is 
represented as 

W = U + T 


where the function U is assumed to be a fixed low order reference field and T 
represents the high order potential. We would expect in this study that dominant 
local gravity anomalies would be specified by a priori values in U and that the 
major problem of fine structure recovery would involve estimating the potential 
T, which would exhibit no local large amplitude anomalous behavior. 


Following report [A] , we choose to represent the high order geopotential by 
local surface density blocks specifying a fictitious surface density layer. 



where G is the gravitational constant, p is the density of the mass layer, and d 
is the distance from the point where T is evaluated to the integration element 
do. The parameter used to describe the high order field is $ = Gp expressed 
in milligals (mgal) , where 

-3 -3 2 

1 mgal = 10 gals = 10 cm/sec 


For order of magnitude comparisons between the surface density parameter $ and 
the expression for corresponding gravity anomalies Ag in milligals, the expression 



is adequate. A more detailed treatment of the surface layer potential is 
given in paper [A] and by Schwarz [15]. 
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For simplicity, the mass layer is considered spread on the surface of a 
sphere with radius equal to the mean radius of the earth. The layer is modeled 
by individual area blocks in which the parameter $ is constant. The disturbing 
potential is then 



j 


til 

where the integration is over the area of the j— block and the sum is over all 

surface blocks. The integral over each block is evaluated by dividing the block 

into a number of sub-blocks of area a.. , where the distance of d.. is the distance 

Jk’ jk 

from the center of the sub-block to the point of evaluation (see Figure 2.1), 

T - 


E 

k 
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The principal portion of this study is concerned with recovery of local 
subsets of blocks near the equator utilizing surface blocks with boundaries 
defined by constant values of latitude and longitude increments (equal-angular 
blocks). The studies involving global recoveries and recoveries near the poles, 
however, utilize equal-area blocks. The algorithm implemented into the GEOMAP 
system for equal-area blocks is principally that due to Morrison [12]. 
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3.0 COMPUTER SOFTWARE 


The GEOMAP/ERODYN software system as developed in report [A] was the basic 
analysis tool utilized for this covariance error analysis study. GEOMAP is a 
large scale orbital parameter estimation and mission analysis program developed 
from the GEODYN program [5], while ERODYN is an orbital and geodetic error 
analysis program designed to be operated directly from GEOMAP output. To 
perform an error analysis of the GEOMAP estimation algorithm, ERODYN requires 
the estimation normal matrix output from GEOMAP. Moreover, to investigate 
by error analysis the partitioning of a parameter set into adjust and qnadjust 
(alias) parameters, the total parameter set must be adjusted in a GEOMAP run 
so that all parameters are included in the normal matrix. This normal matrix, 
obtained from a single GEOMAP run, is then partitioned into appropriate adjust 
and unadjust sets within ERODYN. The setting of a priori parameter sigmas and 
the partitioning of the total parameter set into adjust and unadjust parameters 
within ERODYN permits a great many error analysis runs to be made from a single 
GEOMAP normal matrix. 

The ERODYN program is designed as a single arc error analysis package. The 
multi-arc capability required for the short data arc recovery simulation is 
obtained by concatenating the normal matrices for each arc generated by GEOPMAP. 
For example, let A^ and A^ denote the portion of the normal matrix corresponding 
to the arc dependent parameters for arcs one and two, respectively, and let 
and C 2 denote that portion corresponding to parameters that are common to each 
arc. The two single arc normal matrixes are then 
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and the concatenated matrix for both arcs is 



The large dimensional multi-arc normal matrix is suitable for analysis by the 
linear algebra capability of ERODYN. 

Some modifications to the GEOMAP/ERODYN system were implemented for the 
efficient performance of this study. These modifications included a fully 
general data simulation mode and a multi-satellite SST observability test 
algorithm in GEOMAP and a dynamic partitioning scheme for the normal matrix 
inversion operation in ERODYN to reduce the large computer core requirement 
for the large parameter sets involved in this study. In addition, the single 
arc normal matrix concatenation program was modified to accomodate multi-satellite/ 
multi-arc configurations. Independent verification of the program modification 
was aided by comparing test runs using this software with results obtained from 
an operating BTS program designed to analyze state vector observability obtained 
from simulated SST data [6]. 
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4.0 COVARIANCE ERROR ANALYSIS 


The error analysis for the Bayesian least squares estimator in GEOMAP is 
described in the ERODYN Program Mathematical Description [9] and in report [A]. 
However, for completeness we reproduce the main points of the mathematical 
treatment. 

Let y be an n dimensional vector of observations modeled by the nonlinear 
equation 


y = F (X, S) + e 

where X and S are parameters of the dynamical system. Here X denotes the para- 
neter set to be adjusted (estimated) and S denotes the parameter set which is 
unadjusted, or constrained to constant values in the estimation process but whose 
uncertainties are to be considered in computing the covariance matrix of the 
estimate. The vector e represents the observation noise, which is mean zero and 
statistically independent. Linearizing about a nominal solution of the parameters 
X^ and yields 

y = C 

= A6X + B6S + e 


6X 


6S 


+ e 


where 


6y = y - F (Xn ,S n ) 

«X = X - x N 


6S = S - S 


A = 


N 

9F(X,S) 

ax 


V S N 
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and 


« _ 3F(X,S) 
B " 3S 


^N’ S N 



Here X and S denotes the actual values of the adjusted and unadjusted parameters 
and 6X = X - X represents the difference between the actual and a priori values 

St 3. 

of the adjusted parameters. 


The normal matrix for all parameters is given by 


T 

: wc 


T 

A WA 

T 

A WB 

T 

B WA 

T 

B WB 


where the weight matrix W is the inverse of the observation noise covariance matrix 
and is assumed to be diagonal, i.e. 


-1 T 

W = E(ee ) 



where E denotes expected values and is the ith observation noise sigma. 


The Bayesian least squares estimate of X is given by 


“ T -1 -1 T T -I _1 -1 '■ 

6X = (A WA + P ) A W6y + (A WA + P ) P 6X 

3 3 3 3 


where 
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and 


P = E[(X - X) (X - X) T ] 
a a a 

is the covariance matrix of the a priori estimate X g . Substitutions yield the 
error 


AX = 6X - SX = (A T WA + P 1 )“ 1 (A T WB5S + A T We - p“ 1 6X ) 

3 cl 3. 

showing the three distinct components due to aliasing, measurement noise, and a 
priori uncertainties. 

Under the assumption that 6S, e, and 6X are uncorrelated errors, the covari- 

ct 

ance matrix of the adjusted parameters is given by 

P = E [ (X — X) (X — X) T ] = (A T WA + P -1 )" 1 

a 


+ (a t wa + p 1 ) 1 (a t wb) V (b t wa)(a t wa + p 1 ) 1 

a b a 


where Vg is the diagonal covariance matrix of the unadjusted parameters, 

V = E(6S6S T ) . 

u 


The sensitivity matrix of adjusted to unadjusted parameters is 


S = = (A T WA + P _1 ) VWB 

96b a 


so that the aliasing error to the ith adjusted parameter due to the jth unadjusted 
parameter is 


a . . 
ij 




The ith diagonal element of P, or the variance of the ith adjusted parameter is 
then given by 
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a 2 . = [~(A T WA + P" 1 ): 1 ] +y^(5..a ] 

x L a xx J Z — r xj S. 


so that the standard deviation of the ith adjusted parameter is the root sum square 
(RSS) of the standard deviations due to data noise and a priori uncertainties 
and the standard deviations due to each unadjusted parameter. 
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TABLE 4.1 


Geopause : 


Geostationary: 


GRAYS AT : 


r = 29431200. m 

h = 23,053.045 km. 

e = 0 
i = 90° 


r = 42164541.3 m 

h = 35,786.386 km. 
e = .000388 
i = 1.01° 


r = 6628155.0 m 
h = 250 km. 


e = 0 
i = 90° 
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The data rate assumed in this error analysis study for local recovery is 
one observation every five seconds along the satellite path and four satellite 
passes per surface density block for all cases using the 5° x 5° set of equal- 
angular blocks centered on the equator. Short arc passes are ascending and 
decending. The spacing of satellite passes is slightly different for the 
2 1/2° x 2 1/2° density block and the 5° x 5° equal-area density block studies 
and will be fully described in Section 6.0. Each data arc is treated as being 
independent. The GEOMAP program (report [A]) generates each arc normal matrix 
separately and then for a given set of short arcs the arc normal matrices are 
concatenated into a single normal matrix for error analysis purposes. 

The SST range-rate signal signature for short arc passes of the "low-low” 
GRAVSAT, the GRAV SAT / GEOPAUSE and the GRAVS AT/GEOSTATIONARY as the GRAVSAT 
passes over a 5° x 5°, 1 mgal surface density block located on the equator is 
shown in Figures (4. 1-4. 6). The marker (A) along the time axis denotes the 
time that the low satellite is directly over the center of the density block. 

In these figures all passes are ascending arcs with the low satellite 15° 
below the equator at epoch. The angle in parenthesis designated AM in the legend 
on Figures 4. 2-4. 5 represents the mean anomaly difference of the h;igh GEOPAUSE 
and low GRAVSAT at the epoch, in degrees. Thus, (-40 p ) indicates that at 
epoch the low satellite is 40° ahead and moving away from the high satellite. 

The signal strength is seen to be greater for the "high-low" configuration 
than the "low-low" configuration. This will be displayed in the error analysis 
cases of Section 6.0. Moreover, the magnitude and shape of the "high-low" 
signal varies with the geometry of the high and low orbit at epoch. The data 
arcs generated for the error analysis studies of Section 6.0 were designed to 
be a realistic mix of geometries for the high and low orbits. This will be 
discussed more fully in Section 6.0. 
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FIGURE k . 1 "Low-Low" GRAVSAT SST Signal Over A 1 MGAL 5° x 5° Density Block. 
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6.0 DATA SIMULATION AND ERROR ANALYSIS 


As demonstrated by Hajela [8] and in report [A], if a local set of density 
blocks is to be estimated by SST data the blocks near the edges will be badly 
corrupted (aliased) by nearby unadjusted blocks, whereas blocks near the center 
of the set may be sufficiently removed from the unadjusted blocks that their 
estimates may be acceptable. Thus, to estimate a given set of blocks it will be 
necessary to simultaneously estimate the set of interest plus a surrounding set 
of "buffer" blocks, whose estimates will be discarded due to aliasing. To 
determine the relationships between estimation accuracy and the estimation 
region, data region, data noise, a priori uncertainties in adjusted and 
unadjusted density blocks and a priori uncertainties in the satellite states we 
follow the approach of report [A] and utilize the method of covariance error 
analysis outlined in Section 4.0. It is recalled that in a covariance analysis only 
the covariance matrix for the estimator is calculated - the least squares estimation 
is postulated and not actually performed. Parameters in the "solve for" or adjust 
mode are assumed to be estimated in the postulated least squares adjustment, while the 
parameters in the "consider" or unadjust mode are assumed to affect the functionality 
between the adjusted parameters and the observations but are not estimated. 

The error analysis results are organized into six sets: 


SET I : The "high-low" GRAVS AT /GEOPAUSE configuration simulated 

with twenty-one short data arcs over 225 5° x 5° equal 
angular blocks arranged in a 15 block by 15 block square 
centered on the equator with data coverage over the center 
25 blocks. 


SET II : The "high-low" GRAVSAT/GEOSTATIONARY configuration 

simulated with twenty-one short data arcs over the same 
density block configuration as SET I. The GEOSTATIONARY 
satellite is positioned at the center of the local set of blocks. 


SET III : The "high-high-low" GRAVSAT/GEOSTATIONARY configuration 

simulated with twenty-one short data arcs over the same density 
block configuration as SET I. The two GEOSTATIONARY satellites 
are positioned 60° to the east and west of the center of the 
local set of blocks and take independent, simultaneous SST range-rate 
measurements of the GRAVSAT. 
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SET IV: 


The "high-low" GRAVSAT/ GEOSTATIONARY configuration simulated 
with ten short data arcs over 225 2 1/2° x 2 1/2° equal angular 
blocks arranged in a 15 block by 15 block square centered on 
the equator with data coverage over the center 25 blocks. 

SET V : The "high-high-low" GRAVSAT/ GEOSTATIONARY configuration simulated 

with twelve short data arcs over 185 5° x 5° equal-area blocks 
arranged in 8 latitude bands centered on the north pole with 
data coverage over the first three bands about the pole. The 
two GEOSTATIONARY satellites are positioned 180° apart. 

SET IV : The "low-low" GRAVSAT configuration simulated with twenty-one 

short data arcs over the same density block configuration as SET I. 

The total parameter set for the error analysis consists of six state para- 
meters for each satellite for each arc plus the entire set of surface density 
parameters. The sensitivity matrix C, consisting of the partial derivatives of 
the observations with respect to the parameter set, is evaluated along a 
nominal trajectory where the a priori values of all surface densities have 
been set to zero. The resulting normal matrix is adequate for linear error 
analysis purposes and provides a substantial savings of computer time. 
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6.1 DATA COVERAGE 


The simulated satellite-to-satellite range-rate data for the six sets 
of error analysis studies described in 6.0 all were generated from a two 
body plus J ^ force model. Using the GEOMAP program with an eleventh order 
Cowell predictor-corrector integrator and a 20 second fixed step size to 
integrate the equations of motion. The data, generated at 5 second intervals, 
was then contaminated with random noise of standard deviation .05 cm/sec. The 
5 second data interval produces approximately 15 observations per equatorial 
5° x 5° equal-angular denisty block per arc. However, the high satellite 
geometry, the data arc spacing, and the geographic region covered by data 
varied among the six sets. The particular features of the data arcs for 
each of the sets will be described in this subsection. It is to be noted that 
the simulation is idealized in that there are no transmission delays at the 
satellites, all measurements are instantaneous , and the drag compensation system 
of the GRAVSAT has been assumed error free. This means that error analysis 
results will be somewhat optimistic, but will serve to determine feasibility. 
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SET I: 


The data coverage for the coplanar, circular, polar orbiting GRAVSAT/ GEOPAUSE 
configuration consists of 21 short arc data passes of the low satellite over 
the center 25 blocks of the 225 5° x 5° equal-angular blocks centered on the 
equator. The passes are alternately ascending and decending with approximately 
1.2° spacing, giving approximately four data passes per block. The ground 
track of the low GRAVSAT over the data region is displayed in Figure 6.1. The 
epoch for each data arc is at -15° latitude for ascending and +15° latitude 
for descending passes, while the data region extends between -12 1/2° and +12 1/2° 
latitude. Because the period of the GRAVSAT at 250 km altitude is approximately 
90 minutes and the period of the GEOPAUSE is approximately 14 hours, there is a 
highly variable geometry possible between the high and low satellites. Using 
a line of sight observability test between the satellites, a realistic mix of 
configurations was calculated which would provide the desired data coverage. 

The orientation of the high satellite relative to the low satellite for the 
21 passes is described in Table 6.1. As defined in Section 5.0, AM represents 
the mean anomaly difference of the high GEOPAUSE and low GRAVSAT at the epoch, 
in degrees. A positive AM indicates that the low satellite is "moving toward" 
the high satellite, while a negative AM indicates that the low satellite is 
"moving away" . 

Two sets of 21 data arcs identical to the above coplanar configuration 
with the exception that the orbit plane of the GEOPAUSE was offset from the 
plane of the GRAVSAT by 30° and 45°, respectively, were generated to study 
the effects of adding out-of-plane components in the data. 
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Table 6.1 Orientation of the 21 Data Arcs for the 
GRAVSAT/GEOPAUSE Configuration 

+ means Low satellite moving toward High satellite 
means Low satellite moving away from High satellite 


Arc 


A - Ascending 

D - Descending AM (at epoch) 


1 

2 

3 

4 

5 

6 

7 

8 
9 

10 

11 

12 

13 

14 

15 


16 

17 

18 

19 

20 
21 


A 

D 

A 

D 

A 

D 

A 

D 

A 

D 

A 

D 

A 

D 

A 

D 

A 

D 

A 

D 

A 


-40 
-50 
+60 
+55 
+90 
+16 
-40 
+ 4 
+50 
-40 
+10 
+55 
- 5 
+ 4 
+90 
+16 
+50 
+55 
+10 
+40 
+90 


in degrees 
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SET II: 


The data coverage for the GRAVSAT/GEOSTATIONARY "high-low" configuration 
consists of 21 short arc data passes of the low satellite over the center 25 
blocks of the 225 5° x 5° equal-angular blocks centered on the equator. The 
epochs and ground tracks of the low satellite are identical to SET I, described 
above and displayed in Figure 6.1. The high geostationary satellite is located 
at the center of the 25 block data region. 

SET III: 


The data coverage for the "high-high-low" GRAVSAT/GEOSTATIONARY configuration 
is identical to that of SET II with the exception that there are two simulated 
geostationary satellites located 60° to the east and west of the center of the 
25 blocks covered by data, respectively. All 21 passes of the GRAVSAT are 
simultaneously observable by both high tracking satellites, and SST range-rate 
observations to the low satellite are generated independently. 

SET IV : 

The data coverage for the GRAVSAT/GEOSTATIONARY "high-low" configuration 
consists of 10 short arc data passes of the low satellite over the center 25 
blocks of 225 2 1/2° x 2 1/2° equal-angular blocks centered on the equator. 

The passes are alternately ascending and descending with approximately 1.2° 
spacing, giving approximately two data passes per block. The ground track 
of the low GRAVSAT over the data region is displayed in Figure 6.2. The 
epoch for each data arc is at -15° latitude for ascending and + 15° latitude 
for descending passes, while the data region extends between -6 1/4° and 
+6 1/4° latitude. The high geostationary satellite is located at the center 
of the 25 block data region. 
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LATITUDE 


FIGURE 6.1 Ground Track of 21 GRAVSAT Data Arcs Over 25 5° x 5° Equa 1 -Angular 
Density Blocks. 



LONGITUDE 



LATITUDE 


FIGURE 6.2 Ground Track of 10 GRAVSAT Data Arcs Over 25 2i° x 2$° Equal Angular 
Density Blocks. 
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SET V: 


The data coverage for the GRAVSAT/GEOSTATIONARY "high-high-low" 
configuration consists of 12 short arc data passes of the low satellite over 
the first three latitude bands of 5° x 5° equal-area blocks centered on the north 
pole. The passes are alternately ascending and descending with approximately 
15° spacing, giving 2 data passes per block in the second latitude band. The 
ground track of the low GRAVSAT over the data region is displayed in Figure 6.3. 
The epoch for each arc is at 70° latitude, while the data region extends from 
75° latitude to the pole. The two relay satellites are located 180° apart 
at 90° and 270° longitude, respectively. The data was simulated such that 
the low satellite was observable only to the relay satellite at 90° longitude 
when its longitude was between 0° and 180°, and to the relay satellite at 
270° longitude when its longitude was between 180° and 360°. Note that the 
polar cap receives approximately 6 times more data than the blocks in the 
second latitude band (See Figure 6.3), while the second latitude band receives 
approximately 2 times more data than the blocks in the third latitude band. 

SET VI : 

The data coverage for the GRAVSAT "low-low" configuration (see report [A]) 
consists of 21 short arc data passes of the two low satellites, separated by a 
geometric angle of 6° at 250 km. altitude, over the center 25 blocks of the 225 
5° x 5° equal-angular blocks centered on the equator. All data was generated 
on ascending passes with approximately 1.2° spacing, giving approximately four 
passes per block. The epoch for each data arc is at -15° latitude, while the 
data region extends between -12 1/2° and +12 1/2° in latitude. The data passes 
are symmetrically defined by the ground track of the midpoint between the two 
satellites. 
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90° 



270 ° 


FIGURE 6,3 Ground Track of 12 Data Arcs For GRAVSAT - 2 GEOSTATIONARY Polar 

Configuration of 5° Equal-Area Blocks. The Geostationary Satellites 
Are Located at 90° and 270° Longitude. 
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6.2 ERROR ANALYSIS STUDY 


The error analysis results were generated under the assumption that a 
priori satellite epoch state errors reflect the orbital knowledge obtained 
from other tracking means. No ground tracking data was included in . the study. 
Assumptions concerning the SST range-rate data and the geometry of the satellite 
configurations were discussed in Section 6.0 and 6.1. Any assumptions that are 
particular to any of the six error analysis sets will be described individually 
in the discussion of the error analysis results of those sets. 

In the tables of solution sets I-VI displaying error sigmas, the dashed 
perimeter encloses the data region, while the set of blocks containing two 
numbers are the adjusted density blocks. The two numbers within each adjusted block 
represent the total RSS error sigma for that block and the error sigma for that 
block due to aliasing by unadjusted blocks in milligals (mgal) , i.e. 


°T0TAL / ° ALIAS . 


The error sigma due to data noise and a priori parameter uncertainties is then 


°N0ISE 


TOTAL 


2 

0 ALIAS 


The numbers in unadjusted blocks represent the aliasing contribution of that 

block to the center adjusted block. The signs on these aliasing contributions 

3AX 

reflect the influence of the sensitivity matrix S = . In the tables of 

dAb 

solution sets I-VI displaying correlations, the number within each adjusted 
block is the correlation of that adjusted block with the center adjusted block. 

The results presented in Set VI for the "low-low" configuration are 
extensions of results obtained in report [A] and are included here for purposes 
of comparison with the "high-low" mission. As in the case of the "low-low" 
configuration, it was found necessary to simultaneously solve for the epoch 
state vector of the low satellite along with the density blocks to obtain a 
reasonable recovery for the blocks. The epoch state errors of the high satellite 
produced a much smaller effect on the recovery, however, and studies in Sets I 
and II demonstrate the result of treating the epoch state vector of the high 
satellite as unadjusted. In solution Sets I-VI, however, the high satellite 
state vector will be included with the low satellite state vector in the adjust mode 
unless explicitly stated otherwise. 
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SET I: 


Solution Set I treats the "high-low" GRAV SAT/ GEOPAUSE configuration as 
described in Section 6.1. As in report [A], the center block is the main focus 
of attention in these comparative studies for quantitative aliasing contributions. 
This error analysis set will investigate the effectiveness of reducing aliasing 
errors by adding "buffer layers", the effects of high satellite state errors 
on the solution, and the relation that a priori uncertainty in the adjusted and 
unadjusted density blocks has on the accuracy of the estimation. The coplanar, 

30° and 45° planar offset GRAVSAT/GEOPAUSE configurations are presented in 
Tables I.I-I.XV, Tables I.XVI-I.XIX, and Tables I.XX-I.XXII, respectively. 

Table I. I displays what should be the most favorable solution for the 
coplanar configuration. The data region is buffered with two adjusted layers 
to reduce aliasing effects, the 4 mgal a priori uncertainty on both the adjusted 
and unadjusted blocks represents a reasonable value for 5° x 5 ° density blocks, 
and the satellite state for both the high and low spacecraft are assumed to be 
virtually error free (a priori uncertainties of both position and velocity set 
to the small value .001 for convenience). The resulting aliasing error of the 
center block is small, and the aliasing contribution of unadjusted blocks 
decreases rapidly with distance from the adjusted set. The dominant error 
source is the "noise plus a priori" contribution arising from 

(a t wa + p ^ 1 ) - 1 . 

The reduction of the aliasing effect due to unadjusted blocks by estimating 
"buffer layers" is demonstrated in Tables I.II-I.III. Here the satellite states 
are again assumed virtually error free, but the a priori uncertainty on both 
the adjusted and unadjusted blocks is set to the large value of 50 mgal. 

This greatly increases the aliasing of the adjusted set so that the effective 
strategy of reducing aliasing errors in the central region of interest by 
estimating "buffer layers" not covered by data can be clearly seen. The error 
of 4.55 mgal for the center block of Table I. II is almost entirely due to the 
aliasing contribution, as is the .76 mgal of Table I. Ill when a single layer is 
introduced. However, the double "buffer layer" of Table I. IV results in .33 
mgal error composed of nearly equal "aliasing" and "noise plus a priori" components. 
For realistic values of unadjusted block a priori, as in Table I. I, the "noise 
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plus a priori" component will dominate the total error, showing the two 
"buffer layers" for the 5° x 5° density blocks to be an effective strategy to 
eliminate aliasing errors in the local estimation process. It is to be 
noted that the "high-low" configuration produces a less sensitive range-rate 
data type to this aliasing problem than the "low-low" concept investigated in 
report [A]. In that case two "buffer layers" were required also, (compare 
Tables I. I and VI. I) but the addition of a single layer for the "high-low" 
produced a much better solution with respect to aliasing than did the addition 
of a single layer for the "low-low". 


Tables I.V-I.XI display the recovery errors to be expected utilizing 
realistic values for the uncertainties of density blocks (4 mgal) and the 
epoch states for the high and low spacecraft. Table I.V uses state uncertainty 
values obtained by Koch [11] for the coplanar GRAVSAT/GEOPAUSE mission and finds 
recovery to an accuracy of .2 mgal. To determine the character of the recovery 
as a function of epoch state a priori errors. Tables I.VI-I.XI utilize multiples 
of the set (2m, 5m, 10m) for the position uncertainties in radial, cross track 
and along track (H,C,L) components for the GRAVSAT and the sets (.4m, lm, 2m) 
and (20m, 50m, 100m) for the position uncertainties in HCL coordinates for the 
GEOPAUSE. The velocity uncertainties are calculated from the relations 


a 

v 

r 



v 

CT = — 

v r 
c 


a 

c 



where ( a r » a c » a L ) and (° v >° v ,a v ^ are t * ie position uncertainites and velocity 

r c L 

uncertainties, respectively, in HCL coordinates, and r and v are the corresponding 
magnitudes of the orbital position and velocity vectors. 


Comparing Tables I. IX and I. VI shows that increasing the high satellite epoch 
state uncertainty from (.4m, 2m, lm) to (20m, 100m, 50m) in HCL coordinates 
increases the total error of the center block from .30 mgal to .33 mgal, while 
the aliasing component remains the same. This demonstrates the relative 
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insensitivity of the estimation accuracy to the uncertainty of the high 
satellite orbit (at the 10’s of meters level). Table I. XII treats the high 
satellite epoch state parameters of Table I. IX as unadjusted parameters 
resulting in a total error increase to .39 mgal, while the aliasing error 
component increases from .04 mgal to .26 mgal. The aliasing due to the high 
satellite orbital uncertainty is approximately equal to the "noise plus a 
priori" contribution. 

Tables I. VI (or I. IX) shows that if the low satellite is assumed accurate 
to (2m, 5m, 10m) in HCL coordinates and the high satellite accurate to (.4m, 
lm, 2m) , or (20m, 50m, 100m) then a priori density block errors of 4 mgal 
can be reduced to approximately .3 mgal. If these orbital errors are multiplied 
by a factor of ten, then (Tables I. VII and I.X.) a priori density block errors 
of 4 mgals can be reduced to approximately 1.2 mgal. 

The correlations of the adjusted blocks with the center adjusted block 
for the error analysis studies of Tables I. I, I. IV and I. VI are presented in 
Tables I.XII-I.XV, respectively. These correlations display the resolvability 
of 5° x 5° blocks with the "high- low" configuration for the GRAVSAT at 250 km. 
altitude. 

The offset planes of the GRAVSAT and GEOPAUSE alters the geometry such 
that cross track information is also contained in the SST data. Comparing 
Tables I. XVI and I. XIX for a 30° offset with the corresponding coplanar 
solution of Table I. II and I. VI demonstrates that the additional geometry 
improves the recovery when the a priori uncertainty of the density blocks 
is large (50 mgal) and has virtually no impact when the a priori uncertainties 
are smaller (4 mgal). For example, the error analysis study of Table I. XVII 
was repeated for the coplanar configuration and the total error of the center 
block was found to be .71 mgal, with an aliasing component of .55 mgal. These 
results can be compared with the corresponding GRAVSAT/GEOSTATIONARY case 
of Table II. VIII which falls between the coplanar and 30° planar offset configurations. 

Tables I.XVII-I.XIX demonstrate the effect of changing the magnitude of the 
a priori uncertainty for the unadjusted blocks while keeping the large a priori 
value of 50 mgal for the adjusted blocks, while Tables I.XIX-I.XXI compare 
the properties of "buffer layers" for the 45° planar offset recovery with 
the coplanar case shown in Tables I.II-I.IV. 


- 36 - 






■ 

■ 

■ 

■ 


■ 

■ 

■ 

■ 







■ 

■ 

■ 

■ 

^^2 

■ 

■ 

■ 

■ 







■ 

■ 

■ 


.0003 

■ 

■ 

■ 

■ 







3-98 

.98 

3.91 

T7o¥ 

3.74 

1.21* 

3.77 

1.54 


m 


3.84 

.96 

3.96 

.74 







3.92 

.97 

3.38 

.57 

m 

1 .09 
.16 

1.11 

.21 

1.07 
. 18 

n 

2.88 

.58 

3.88 

.71 

i 

j 






3.85 

Ka 

14.08 ' 

.26 

.26 

■ .27 

EH 

1.62 

3.82 




1 .09 

m 

E9 

.04 

.04 

.04 

.20 

.7^ 




3.78 

TTo^T 



.18 

.05 

.18 

.04 

.19 
. 04 

p| 

El 

3.79 

.70 




OOA 

-005 

-.006 



H 

.17 

.04 

.17 

.04 

.18 

.04 

.29 ■ 
.06 , 

sgsj 

■B 

3.75 

~ 7£5 

-.007 

-.005 





3 . 66 
1.11 

.85 

1.36 

El 

.19 

.05 

.18 

.04 

mm 

• 30 | 

.06 ; 

1.02 

.16 








3-75 

1.18 

.93 

.20 


.29 

.04 

• 27 
.04 

.29 

.06 

2jj| 

wm 

wm 

3.82 
• 76 







3.86 

.95 

2.65 

.55 

■mm 

1.30 

.20 

1 .09 
.15 

1.25 

.22 

1.43 

.36 


3-87 

.85 



j 




3.96 

.85 

3.82 

.98 

3.71 

1.17 

3.66 

1.52 

3-78 

1.72 

3.61 

1.40 

3.63 
1 .20 

3.82 

1.11 

3.95 

.84 







■ 

■ 

■ 

■ 

122 

■ 

■ 

■ 

■ 







■ 

■ 


■ 

8 

■ 

■ 

■ 

■ 







■ 

■ 


■ 


■ 

■ 

■ 

■ 





TABLE l.l GRAVSAT - GEOPAUSE (COPLANAR) 

A POSTERIORI DENSITY SIGMAS 
21 SHORT DATA ARCS (DATA NOISE = .05 CM/SEC) 

ADJUSTED PARAMETERS AND A PRIORI SIGMAS: UNADJUSTED PARAMETERS AND UNCERTAINTIES 
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TABLE l.ll GRAVSAT - GEOPAUSE (COPLANAR) 

A POSTERIORI DENSITY SIGMAS 
21 SHORT DATA ARCS (DATA NOISE = .05 CM/SEC) 

ADJUSTED PARAMETERS AND A PRIORI SIGMAS: UNADJUSTED PARAMETERS AND UNCERTAINTIES 
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TABLE I . I I I GRAVSAT - GEOPAUSE (COPLANAR) 

A POSTERIORI DENSITY SIGMAS 
21 SHORT DATA ARCS (DATA NOISE = .05 CM/SEC) 

ADJUSTED PARAMETERS AND A PRIORI SIGMAS: UNADJUSTED PARAMETERS AND UNCERTAINTIES 
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TABLE I.IV GRAVSAT - GEOPAUSE (COPLANAR) 

A POSTERIORI DENSITY SIGMAS 
21 SHORT DATA ARCS (DATA NOISE = .05 CM/SEC) 

ADJUSTED PARAMETERS AND A PRIORI SIGMAS: UNADJUSTED PARAMETERS AND UNCERTAINTIES 

81 DENSITY BLOCKS, a = 50 MGAL 1M DENSITY BLOCKS, a = 50 MGAL 
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TABLE l.v GRAVSAT - GEOPAUSE (COPLANAR) 

A POSTERIORI DENSITY SIGMAS 
21 SHORT DATA ARCS (DATA NOISE = .05CM/SEC) 

ADJUSTED PARAMETERS AND A PRIORI SIGMAS: UNADJUSTED PARAMETERS AND UNCERTAINTIES 

81 DENSITY BLOCKS, a = 4 MGAL 144 DENSITY BLOCKS, 0=4 MGAL 
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TABLE I. VI GRAVSAT - GEOPAUSE (COPLANAR) 

A POSTERIORI DENSITY SIGMAS 
21 SHORT DATA ARCS (DATA NOISE = .05 CM/SEC) 

ADJUSTED PARAMETERS AND A PRIORI SIGMAS: UNADJUSTED PARAMETERS AND UNCERTAINTIES 

81 DENSITY BLOCKS, a = 4 MGAL 144 DENSITY BLOCKS, a = 4 MGAL 

LOW SATELLITE STATE, a [HCL]= (2M,5M,10M) 

HIGH SATELLITE STATE, a [HCL]= (.4M,IM,2M) 
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TABLE I .VI I GRAVSAT - GEOPAUSE (COPLANAR) 

A POSTERIORI DENSITY SIGMAS 
21 SHORT DATA ARCS (DATA NOISE = .05 CM/SEC) 

ADJUSTED PARAMETERS AND A PRIORI SIGMAS: UNADJUSTED PARAMETERS AND UNCERTAINTIES 

81 DENSITY BLOCKS, a = b MGAL DENSITY BLOCKS, o = k MGAL 

LOW SATELLITE STATE, a [HCL]= (10M,25M,50M) 

HIGH SATELLITE STATE, o [HCL]= (2M,5M,10M) 
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TABLE I .VI 1 1 GRAVSAT - GEOPAUSE (COPLANAR) 

A POSTERIORI DENSITY SIGMAS 
21 SHORT DATA ARCS (DATA NOISE = .05 CM/SEC) 

ADJUSTED PARAMETERS AND A PRIORI SIGMAS: UNADJUSTED PARAMETERS AND UNCERTAINTIES 

81 DENSITY BLOCKS, a = A MGAL 144 DENSITY BLOCKS, o = 4 MGAL 

LOW SATELLITE STATE, a [HCL]= (20M, 50M, 1 OOM) 

HIGH SATELLITE STATE, o [HCL]= (4M,10M,20M) 
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TABLE MX GRAVSAT - GEOPAUSE (COPLANAR) 

A POSTERIORI DENSITY SIGMAS 
21 SHORT DATA ARCS (DATA NOISE = .05 CM/SEC) 

ADJUSTED PARAMETERS AND A PRIORI SIGMAS: UNADJUSTED PARAMETERS AND UNCERTAINTIES 

81 DENSITY BLOCKS, a = 4 MGAL 144 DENSITY BLOCKS, a = 4 MGAL 

LOW SATELLITE STATE, 0 [HCL]= (2M.5M.10M) 

HIGH SATELLITE STATE, o [HCL]= (20M.50M, 100M) 
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TABLE I.X GRAVSAT - GEOPAUSE (COPLANAR) 

A POSTERIORI DENSITY SIGMAS 
21 SHORT DATA ARCS (DATA NOISE = .05 CM/SEC) 

ADJUSTED PARAMETERS AND A PRIORI SIGMAS: UNADJUSTED PARAMETERS AND UNCERTAINTIES 

81 DENSITY BLOCKS, o = 4 MGAL 144 DENSITY BLOCKS, 0=4 MGAL 

LOW SATELLITE STATE, 0 [HCL]= (20M.50M, 100M) 

HIGH SATELLITE STATE, 0 [HCL]= (200M,500M, 1000M) 
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TABLE I .XI GRAVSAT - GEOPAUSE (COPLANAR) 

A POSTERIORI DENSITY SIGMAS 
21 SHORT DATA ARCS (DATA NOISE = .05 CM/SEC) 

ADJUSTED PARAMETERS AND A PRIORI SIGMAS: UNADJUSTED PARAMETERS AND UNCERTAINTIES 

81 DENSITY BLOCKS, a = A MGAL lM DENSITY BLOCKS, a = b MGAL 

LOW SATELLITE STATE, a [HCL]= (40M, 100M.200M) 

HIGH SATELLITE STATE, a [HCL]= (400M, 1000M,2000M) 
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TABLE I .XI I GRAVSAT - GEOPAUSE (COPLANAR) 

A POSTERIORI DENSITY SIGMAS 
21 SHORT DATA ARCS (DATA NOISE = .05 CM/SEC) 

ADJUSTED PARAMETERS AND A PRIORI SIGMAS: UNADJUSTED PARAMETERS AND UNCERTAINTIES 

81 DENSITY BLOCKS, 0=4 MGAL 144 DENSITY BLOCKS, a = 4 MGAL 

LOW SATELLITE STATE, a [HCL]= (2M,5M,10M) 

HIGH SATELLITE STATE, a [HCL]= (20M, 50M, 1 OOM) 
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TABLE I.XIII GRAVSAT - GEOPAUSE (COPLANAR) 

CORRELATIONS OF ADJUSTED BLOCKS WITH CENTER BLOCK 
21 SHORT DATA ARCS (DATA NOISE = .05 CM. SEC) 

ADJUSTED PARAMETERS AND A PRIOR! SIGMAS: UNADJUSTED PARAMETERS AND UNCERTAINTIES 

81 DENSITY BLOCKS, a - k MGAL DENSITY BLOCKS, a = k MGAL 

LOW SATELLITE STATE, o [HCL]= (.001M, .OOIM, .001M) 

HIGH SATELLITE STATE, a [HCL]= ( .OOIM, .OOIM, .OOIM) 
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TABLE I.XIV GRAVSAT - GEOPAUSE (COPLANAR) 

CORRELATIONS OF ADJUSTED BLOCKS WITH CENTER BLOCK 
21 SHORT DATA ARCS (DATA NOISE = .05 CM/SEC) 

ADJUSTED PARAMETERS AND A PRIORI SIGMAS: UNADJUSTED PARAMETERS AND UNCERTAINTIES: 

8) DENSITY BLOCKS, a = 50 MGAL ]kk DENSITY BLOCKS, a = 50 MGAL 

LOW SATELLITE STATE, a [HCL]= ( .001 M, . 001M, . 00 1M) 

HIGH SATELLITE STATE, o [HCL]= ( . 00 1 M, . 001 M, . 001 M) 
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TABLE I.XV GRAVSAT - GEOPAUSE (COPLANAR) 

CORRELATIONS OF ADJUSTED BLOCKS WITH CENTER BLOCK 
21 SHORT DATA ARCS (DATA NOISE = .05 CM/SEC) 

ADJUSTED PARAMETERS AND A PRIORI SIGMAS: UNADJUSTED PARAMETERS AND UNCERTAINTIES 

81 DENSITY BLOCKS, o = 4 MGAL 144 DENSITY BLOCKS, a = 4 MGAL 

LOW SATELLITE STATE,, a [HCL]= (2M,5M,10M) 

HIGH SATELLITE STATE, a [HCL]= (.4M,1M,2M) 
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TABLE I.XVI GRAVSAT - GEOPAUSE (30° Off-Plane) 

A POSTERIORI DENSITY SIGMAS 
21 SHORT DATA ARCS (DATA NOISE = .05 CM/SEC) 

ADJUSTED PARAMETERS AND A PRIORI SIGMAS: UNADJUSTED PARAMETERS AND UNCERTAINTIES 

81 DENSITY BLOCKS, a » 50 MGAL 1M DENSITY BLOCKS, a = 50 MGAL 

LOW SATELLITE STATE, o [HCL]= ( . 001 M, . 00 1 M, . 00 1 M) 

HIGH SATELLITE STATE, o [HCL)= (.001M, .OOTM, .001M) 
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TABLE I . XV I I 


GRAVSAT - GEOPAUSE (30° Off-Plane) 

A POSTERIORI DENSITY SIGMAS 
21 SHORT DATA ARCS (DATA NOISE = .05 CM/SEC) 


ADJUSTED PARAMETERS AND A PRIORI SIGMAS: 

81 DENSITY BLOCKS, a = 50 MGAL 

LOW SATELLITE STATE, a [HCL]= (2M.5M, I0M) 

HIGH SATELLITE STATE, a [HCL]= (.AM.1M.2M) 


UNADJUSTED PARAMETERS AND UNCERTAINTIES 
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TABLE I. XVI II GRAVSAT - GEOPAUSE (30° Off-Plane) 

A POSTERIORI DENSITY SIGMAS 
21 SHORT DATA ARCS (DATA NOISE = .05 CM/SEC) 

ADJUSTED PARAMETERS AND A PRIORI SIGMAS: UNADJUSTED PARAMETERS AND UNCERTAINTIES 

8 1 DENSITY BLOCKS a = 50 MGAL ]kh DENSITY BLOCKS, a = k MGAL 

LOW SATELLITE STATE, a [HCL]= (2M,5M,10M) 

HIGH SATELLITE STATE, a [HCL]= (.4M,1M,2M) 
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TABLE I, XIX GRAVSAT - GEOPAUSE (30° Off-Plane) 

A POSTERIORI DENSITY SIGMAS 
21 SHORT DATA ARCS (DATA NOISE = .05 CM/SEC) 

ADJUSTED PARAMETERS AND A PRIORI SIGMAS: UNADJUSTED PARAMETERS AND UNCERTAINTIES 

81 DENSITY BLOCKS, o = 4 MGAL 1M DENSITY BLOCKS, a = 4 MGAL 

LOW SATELLITE STATE, a [HCL]= (2M,5M,10M) 

HIGH SATELLITE STATE, o [HCL]= (.4M,1M,2M) 
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TABLE I.XX GRAVSAT - GEOPAUSE (^5° Off-Plane) 

A POSTERIORI DENSITY SIGMAS 
21 SHORT DATA ARCS (DATA NOISE = .05 CM/SEC) 

ADJUSTED PARAMETERS AND A PRIORI SIGMAS: UNADJUSTED PARAMETERS AND UNCERTAINTIES 

25 DENSITY BLOCKS, a = 50 MGAL 200 DENSITY BLOCKS, a = 50 MGAL 

LOW SATELLITE STATE, a [HCL]= ( . 001M, .00 1M, .001M) 

HIGH SATELLITE STATE, a [HCL]= ( .001M, . 001M, . 001M) 
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TABLE I, XXI GRAVSAT - GEOPAUSE (45° Off-Plane) 

A POSTERIORI DENSITY SIGMAS 
21 SHORT DATA ARCS (DATA NOISE = .05 CM/SEC) 

ADJUSTED PARAMETERS AND A PRIORI SIGMAS: UNADJUSTED PARAMETERS AND UNCERTAINTIES 

49 DENSITY BLOCKS, a = 50 MGAL 176 DENSITY BLOCKS, a = 50 MGAL 

LOW SATELLITE STATE, a [HCL]= (.001M, .001M, .001M) 

HIGH SATELLITE STATE, a [HCL]= (.001M, .001M,.001M) 
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TABLE I . XX I I GRAVSAT - GEOPAUSE (45° Off-Plane) 

A POSTERIORI DENSITY SIGMAS 
21 SHORT DATA ARCS (DATA NOISE = .05 CM/SEC) 

ADJUSTED PARAMETERS AND A PRIORI SIGMAS: UNADJUSTED PARAMETERS AND UNCERTAINTIES 

81 DENSITY BLOCKS, a = 50 MGAL 144 DENSITY BLOCKS, a = 50 MGAL 

LOW SATELLITE STATE, a [HCL]= ( 1 001 M, . 001M, . OOIM) 

HIGH SATELLITE STATE, a [HCL]= (. 001 M, . OOIM, . OOIM) 















































































SET II: 


Solution Set II investigates the "high-low" GRAVSAT/GEOSTATIONARY 
configuration as described in Section 6.1. This study is designed to 
compare the efficiency of the geostationary high satellite concept for 
gravity find structure recovery with the "low-low" configuration of report [A] 
(Set VI) and the high GEOPAUSE satellite of Set I. 

Table II. I displays the most favorable case for the high GEOSTATIONARY 
with two "buffer layers", a realistic 4 mgal surface density uncertainty and 
nearly perfect knowledge of the epoch states, while Table II. II shows the same 
situation with 50 mgal uncertainty for the density blocks. Comparison with 
solution Set I shows that the GEOPAUSE and GEOSTATIONARY spacecraft have 
very similar characteristics with regard to recovery of gravity fine structure. 
The geostationary configuration has a slight advantage when the density block 
a priori uncertainties are small (4 mgal) and the GEOPAUSE configuration 
is more efficient when the a priori uncertainties are very large (50 mgal) . 

The case displayed in Table II. I, assuming error free satellite epoch 
state errors can be compared with results obtained by Hajela [8] for a "high- 
low" configuration with the low satellite at 250 kn. altitude. His solution 
set 250-5-2 for the estimation of 24 5° x 5° equal-area gravity anomaly blocks, 
with data over the interior 12 blocks, yielded a R.M. S. standard deviation 
of 3.4 mgal for the interior anomaly blocks. Hajela assumed a data noise 
of .08 cm/sec, a data interval of 10 seconds, and two passes per block. His 
result can be compared to the result of the center density block of Table II. I 
by the adjustment 

(^)-( 3 .4 ■ .17 mgal 

This compares very closely with the .15 mgal found in Table II. I. 

As in Set I, to determine the character of the density block recovery as 
a function of epoch state a priori errors, Tables II.III-II.VII utilize multiples 
of (2m, 5m, 10m) and (20m, 50m, 100m) in HCL coordinates for the GRAVSAT and 
GEOSTATIONARY position uncertainties, respectively. These uncertainties, with a 
4 mgal uncertainty for the density blocks results in an error of .27 mgal 
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in the recovery of the center block, while increasing the state a priori 
by a factor of 20 results in an error of 1.36 mgal. 

Tables II. VIII and II. IX, together with II. IV demonstrate the effect of 
changing the magnitude of the a priori uncertainty for the unadjusted blocks 
from 50 mgal to 4 mgal while keeping the large a priori value of 50 mgal for 
the adjusted blocks. This is to be compared with Tables I.XIII-I.XV. 

The effect of high satellite unadjusted state errors contaminating the 
recovery is observed in Table II. X. Here geostationary state errors of 
(20m, 50m, 100m) and 4 mgal uncertainty in these density blocks result in 
an error in the center adjusted block of .29 mgal, with an aliasing component 
of .14 mgal. This is a slightly degraded solution relative to that of Table II. IV. 
The recovery of gravitational fine structure from SST range-rate is seen to be 
insensitive to the errors of the high satellite when these errors are reasonably 
small (order of tens of meters) . 

The effect of reduced data noise is displayed in Table II. XI where the 
noise sigma is set to .01 cm/sec, a factor of fine improvement, resulting in 
a total recovery error of .22 mgal. This is to be compared with the solution 
of Table II. IV with a noise sigma of .05 cm/sec, where we found a recovery 
error of .27 mgal. The aliasing component of error is the same in both instances, 
showing the dominating error source to be the a priori uncertainties in the 
adjusted parameters. 

The GRAVSAT satellite has been presented in these studies as a drag free 
satellite with a surface force compensation system and thrusters to maintain 
an orbit due to only gravitational forces. However, following the error analysis 
case study of Table II. V, drag and drag-rate parameters were added to the set 
of unadjusted parameters for the low satellite for each of the twenty-one arcs. 

The error analysis displayed virtually no impact of the drag parameters (even 
with large a priori errors) on the recovery of the density blocks. This result 
is strongly dependent, of course, on the assumed use of short (approximately 
10 minutes) , independent data arcs with reasonably small a priori uncertainties 
of the low satellite epoch state vector. 
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The correlations of the adjusted blocks with the center block for the error 
analysis studies of Tables II. I and II. IV are presented in Tables II. XIII and 
II. XIV, respectively. The correlations display good resolvability in the recovery 
of the 5 ° x 5° blocks, and compare very closely with the correlations for the 
corresponding error analysis studies of the GRAVSAT/GEOPAUSE configuration in 
Set I. 
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TABLE I I . I 


GRAVSAT - GEOSTATIONARY 

A POSTERIORI DENSITY SIGMAS 
21 SHORT DATA ARCS (DATA NOISE = .05 CM/SEC) 


ADJUSTED PARAMETERS AND A PRIORI SIGMAS: 

81 DENSITY BLOCKS, a = k MGAL 

LOW SATELLITE STATE, a [HCL]= ( . 001 M, . 001 M, . 001 M) 
HIGH SATELLITE STATE, a [HCL]= ( . 001 M, . 001 M, . 001 M) 


UNADJUSTED PARAMETERS AND UNCERTAINTIES 
]kk DENSITY BLOCKS, = 4 MGAL 
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TABLE 11.11 GRAVSAT - GEOSTATIONARY 

A POSTERIORI DENSITY SIGMAS 
21 SHORT DATA ARCS (DATA NOISE = .05 CM/SEC) 

ADJUSTED PARAMETERS AND A PRIORI SIGMAS: UNADJUSTED PARAMETERS AND UNCERTAINTIES 

8 1 DENSITY BLOCKS, a = 50 MGAL DENSITY BLOCKS, o = 50 MGAL 

LOW SATELLITE STATE, a [HCL]= ( . 001 M, . 001 M, . 001 M) 

HIGH SATELLITE STATE, a [HCL]= ( . OOIM, . 001M , . 001M) 
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TABLE I I . ! I I 


GRAVSAT - GEOSTATIONARY 

A POSTERIORI DENSITY SIGMAS 
21 SHORT DATA ARCS (DATA NOISE = .05 CM/SEC) 


ADJUSTED PARAMETERS AND A PRIORI SIGMAS: 

81 DENSITY BLOCKS, a = A MGAL 

LOW SATELLITE STATE, o [HCL]= (1M,2.5M,5M) 

HIGH SATELLITE STATE, a [HCL]= ( 1 0M.25M.50M) 


UNADJUSTED PARAMETERS AND UNCERTAI NTI ES 
DENSITY BLOCKS,, a = k MGAL 
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TABLE II. IV GRAVSAT - GEOSTATIONARY 

A POSTERIORI DENSITY SIGMAS 
21 SHORT DATA ARCS (DATA NOISE = .05 CM/SEC) 

ADJUSTED PARAMETERS AND A PRIORI SIGMAS: UNADJUSTED PARAMETERS AND UNCERTAINTIES 

81 DENSITY BLOCKS, a = 4 MGAL 144 DENSITY BLOCKS, a = 4 MGAL 

LOW SATELLITE STATE, a [HCL]= (2M.5M.10M) 

HIGH SATELLITE STATE, a [HCLJ= (20M, 50M, IOOM) 
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TABLE II.V GRAVSAT - GEOSTATIONARY 

A POSTERIORI DENSITY SIGMAS 
21 SHORT DATA ARCS (DATA NOISE = .05 CM/SEC) 

ADJUSTED PARAMETERS AND A PRIORI SIGMAS: UNADJUSTED PARAMETERS AND UNCERTAINTIES: 

81 DENSITY BLOCKS, a = 4 MGAL 144 DENSITY BLOCKS, cr = 4 MGAL 

LOW SATELLITE STATE, o [HCL]= (10M,25M,50M) 

HIGH SATELLITE STATE, a [HCL]= (100M,250M,500M) 
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TABLE II. VI GRAVSAT - GEOSTATIONARY 

A POSTERIORI DENSITY SIGMAS 
21 SHORT DATA ARCS (DATA NOISE = .05 CM/SEC) 

ADJUSTED PARAMETERS AND A PRIORI SIGMAS: UNADJUSTED PARAMETERS AND UNCERTAINTIES 

81 DENSITY BLOCKS, a = k MGAL 144 DENSITY BLOCKS, a = k MGAL 

LOW SATELLITE STATE, a [HCL]= (20M,50M, 1 OOM) 

HIGH SATELLITE STATE, a [HCL]= (200M,500M, 1000M) 
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TABLE II. VI I GRAVSAT - GEOSTATIONARY 

A POSTERIORI DENSITY SIGMAS 
21 SHORT DATA ARCS (DATA NOISE = .05 CM/SEC) 

ADJUSTED PARAMETERS AND A PRIORI SIGMAS: UNADJUSTED PARAMETERS AND UNCERTAINTIES 

81 DENSITY BLOCKS, a = k MGAL ]kk DENSITY BLOCKS, a = b MGAL 

LOW SATELLITE STATE, a [HCL]= (40M, 100M,200M) 

HIGH SATELLITE STATE, a [HCL]= (AOOM, 10C0M.2000M) 
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TABLE II. VI I I GRAVSAT - GEOSTATIONARY 

A POSTERIORI DENSITY SIGMAS 
21 SHORT DATA ARCS (DATA NOISE = .05 CM/SEC) 

ADJUSTED PARAMETERS AND A PRIORI SIGMAS: UNADJUSTED PARAMETERS AND UNCERTAINTIES 

81 DENSITY BLOCKS, a = 50 MGAL 14A DENSITY BLOCKS, a = 50 MGAL 

LOW SATELLITE STATE, a [HCL]= (2M,5M,10M) 

HIGH SATELLITE STATE, a [HCL]= (20M.50M, 100M) 
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TABLE II. IX GRAVSAT - GEOSTATIONARY 

A POSTERIORI DENSITY SIGMAS 
21 SHORT DATA ARCS (DATA NOISE = .05 CM/SEC) 

ADJUSTED PARAMETERS AND A PRIORI SIGMAS: UNADJUSTED PARAMETERS AND UNCERTAINTIES 

81 DENSITY BLOCKS, a = 50 MGAL 144 DENSITY BLOCKS, a = k MGAL 

LOW SATELLITE STATE, a [HCL]= (2M,5M,10M) 

HIGH SATELLITE STATE, a [HCL]= (20M.50M.I00M). 
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GRAVSAT - GEOSTATIONARY 

A POSTERIORI DENSITY SIGMAS 
21 SHORT DATA ARCS (DATA NOISE = .05 CM/SEC) 


ADJUSTED PARAMETERS AND A PRIORI SIGMAS: 

81 DENSITY BLOCKS, o = k MGAL 

LOW SATELLITE STATE, a [HCL]= (2M.5MJ0M) 
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TABLE X.XI GRAVSAT - GEOSTATIONARY 

A POSTERIORI DENSITY SIGMAS 
21 SHORT DATA ARCS (DATA NOISE = .01 CM/SEC) 

ADJUSTED PARAMETERS AND A PRIORI SIGMAS: UNADJUSTED PARAMETERS AND UNCERTAINTIES 

81 DENSITY BLOCKS, o = 4 MGAL 144 DENSITY BLOCKS, a = 4 MGAL 

LOW SATELLITE STATE, a [HCL]= (2M.5M.10M) 

HIGH SATELLITE STATE, 0 [HCL]= (20M,50M, JOOM) 
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TABLE II. XI I GRAVSAT - GEOSTATIONARY (Drag S Drag-Rate Unadjusted) 

A POSTERIORI DENSITY SIGMAS 
21 SHORT DATA ARCS (DATA NOISE => .05 CM/SEC) 

ADJUSTED PARAMETERS AND A PRIORI SIGMAS: UNADJUSTED PARAMETERS AND UNCERTAINTIES 

81 DENSITY BLOCKS, a = 4 MGAL DENSITY BLOCKS, a = U MGAL 

LOW SATELLITE STATE, 0 [HCL]= (10M.25M.50M) 

HIGH SATELLITE STATE, a [HCL]= (100M,250M,500M) 
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TABLE II.XIII GRAVSAT - GEOSTATIONARY 

CORRELATIONS OF ADJUSTED BLOCKS WITH CENTER BLOCK 
21 SHORT DATA ARCS (DATA NOISE - .05 CM/SEC) 

ADJUSTED PARAMETERS AND A PRIORI SIGMAS: UNADJUSTED PARAMETERS AND UNCERTAINTIES 

81 DENSITY BLOCKS, a = A MGAL 1 hh DENSITY BLOCKS, o = b MGAL 

LOW SATELLITE STATE, a [HCL]= (.001M..001M, .001M) 

HIGH SATELLITE STATE, o [HCL]= (. 001 M, .001 M, . 001 M) 
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TABLE II. XIV GRAVSAT - GEOSTATIONARY 

CORRELATIONS OF ADJUSTED BLOCKS WITH CENTER BLOCK 
21 SHORT DATA ARCS (DATA NOISE = .05 CM/SEC) 

ADJUSTED PARAMETERS AND A PRIORI SIGMAS: UNADJUSTED PARAMETERS AND UNCERTAINTIES 

81 DENSITY BLOCKS, o = 4 MGAL 144 DENSITY BLOCKS, a = 4 MGAL 

LOW SATELLITE STATE, a [HCL]= (2M.5M.10M) 

HIGH SATELLITE STATE, a [HCL]= (20M.50M, 100M) 
















































SET III: 


Solution Set III investigates the "high-high-low" GRAVSAT/GEOSTATIONARY 
configuration utilizing two high spacecraft described in Section 6.1. The 
two high satellites, 60° to the east and west of the center block, simultaneously 
observe the low satellite throughout all twenty-one arcs, introducing more cross 
track geometry in the data type and resulting in twice the data coverage of the 
single high geostationary satellite of solution Set II. Thus, in comparing 
with Sets I and II, the error analysis results in Tables III.I-III.VI should be 
multiplied by 

Tables III.I-III.VI repeat the studies of Tables II.II-II.VII for a single 
high spacecraft over the center of the data region. Considering the factor, 
we find the same result observed when comparing the coplanar and planar offset 
GRAVS AT/ GEOPAUSE configurations of Set I; for large values of the a priori uncer- 
tainty of unadjusted density blocks (50 mgal) the "out-of-plane" geometry of the 
two high satellites (or planar offset GEOPAUSE) provides a better recovery, while 
for a smaller more reasonable value of unadjusted density block a priori for fine 
structure estimation (4 mgal) the more "in plane" configuration provides the 
better recovery. This is due to the fact that for short arcs, although the low 
satellite is most strongly perturbed by the unadjusted blocks nearest the 
ground track, large uncertainties in unadjusted blocks away from the ground 
track represent what can be substantial "out-of-plane" orbital perturbations, 
while for small uncertainties of the unadjusted blocks, the orbital perturbations 
are almost totally planar. 

To illustrate this effect. Tables III.VII-III.XII repeat the studies of 
Tables III.I-III.VI, but with only the single high satellite 60° to the west of 
the center block observing the GRAVSAT. The total data used is then half that 
of Tables III.I-III.VI, and exactly the same as the "high-low" studies of 
Sets I and II. Comparing with the "high-low" GRAVSAT/GEOSTATIONARY case 
of Set II (with 4 mgal a priori uncertainties on the density blocks and the 
high satellite directly over the center block) we find that, for the smaller 
epoch state errors, the recovery uncertainties of Tables III.VII-III.X are nearly 
twice those of Set II. The recovery uncertainties for the large epoch state 
errors are similar in magnitude. Multiplying the uncertainties of Tables III.I- 
III.VI by \ fZ and comparing with the uncertainties of Tables III.VII-III.XII 
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shows that the added geometry of the second high satellite 60° east of center 
improves the solution uncertainty of the single high satellite 60° west of 
center by approximately 20%. 
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TABLE ill, I GRAVSAT - 2 GEOSTATIONARY 

A POSTERIORI DENSITY SIGMAS 
21 SHORT DATA ARCS (DATA NOISE = .05 CM/SEC) 

ADJUSTED PARAMETERS AND A PRIORI SIGMAS: UNADJUSTED PARAMETERS AND UNCERTAINTIES 

8] DENSITY BLOCKS, o = 50 MGAL ]kk DENSITY BLOCKS, a = 50 MGAL 

LOW SATELLITE STATE, a [HCL]= ( . 001 M , . 001 M, .001M) 

HIGH SATELLITE STATE, a [HCL]= (.OOIM, .OOIM, .OOIM) 
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TABLE I I I . I I GRAVSAT - 2 GEOSTATIONARY 

A POSTERIORI DENSITY SIGMAS 
21 SHORT DATA ARCS (DATA NOISE = .05 CM/SEC) 

ADJUSTED PARAMETERS AND A PRIORI SIGMAS: UNADJUSTED PARAMETERS AND UNCERTAINTIES 

81 DENSITY BLOCKS, a = b MGAL 144 DENSITY BLOCKS, a = 4 MGAL 

LOW SATELLITE STATE, a [HCL]= (1M,2.5M,5M) 

HIGH SATELLITE STATE, a [HCL]= (10M,25M,50M) 
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TABLE I I 1,111 GRAVSAT - 2 GEOSTATIONARY 

A POSTERIORI DENSITY SIGMAS 
21 SHORT DATA ARCS (DATA NOISE = .05 CM/SEC) 

ADJUSTED PARAMETERS AND A PRIORI SIGMAS: UNADJUSTED PARAMETERS AND UNCERTAINTIES 

81 DENSITY BLOCKS, a = b MGAL 1 DENSITY BLOCKS, a = b MGAL 

LOW SATELLITE STATE, a [HCL]= (2M.5M.10M) 

HIGH SATELLITE STATE, a [HCL]= (20M.50M, 100M) 
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TABLE I I!. IV GRAVSAT - 2 GEOSTATIONARY 

A POSTERIORI DENSITY SIGMAS 
21 SHORT DATA ARCS (DATA NOISE = .05 CM/SEC) 

ADJUSTED PARAMETERS AND A PRIORI SIGMAS: UNADJUSTED PARAMETERS AND UNCERTAINTIES 

81 DENSITY BLOCKS, a = 4 MGAL 144 DENSITY BLOCKS, a = 4 MGAL 

LOW SATELLITE STATE, a [HCL]= ( 1 OM, 25M.50M) 

HIGH SATELLITE STATE, a [HCL]= (!00M,250M,500M) 
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TABLE lll.v GRAVSAT - 2 GEOSTATIONARY 

A POSTERIOR! DENSITY SIGMAS 
21 SHORT DATA ARCS (DATA NOISE = .05 CM/SEC) 

ADJUSTED PARAMETERS AND A PRIORI SIGMAS: UNADJUSTED PARAMETERS AND UNCERTAINTIES 

81 DENSITY BLOCKS, a = 4 MGAL 144 DENSITY BLOCKS, a = 4 MGAL 

LOW SATELLITE STATE, a [HCL]= (20M.50M, 100M) 

HIGH SATELLITE STATE, a [HCL]= (200M.500M, 1000M) 
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TABLE I II. VI GRAVSAT - 2 GEOSTATIONARY 

A POSTERIORI DENSITY SIGMAS 
21 SHORT DATA ARCS (DATA NOISE = ,05 CM/SEC) 

ADJUSTED PARAMETERS AND A PRIORI SIGMAS: UNADJUSTED PARAMETERS AND UNCERTAINTIES 

8! DENSITY BLOCKS, a = b MGAL DENSITY BLOCKS, a = k MGAL 

LOW SATELLITE STATE, a [HCL]= (^OM, I00M.200M) 

HIGH SATELLITE STATE, a [HCL]= (400M, 1 000M,2000M) 
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TABLE I II. VIII GRAVSAT - I GEOSTATIONARY (60° West of Center) 

A POSTERIORI DENSITY SIGMAS 
21 SHORT DATA ARCS (DATA NOISE = .05 CM/SEC) 

ADJUSTED PARAMETERS AND A PRIORI SIGMAS; UNADJUSTED PARAMETERS AND UNCERTAINTIES 

8 1 DENSITY BLOCKS, o = k MGAL \kk DENSITY BLOCKS, a = k MGAL 

LOW SATELLITE STATE, a [HCL]= (1M,2.5M,5M) 

HIGH SATELLITE STATE, a [HCL]= (10M,25M,50M) 
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TABLE I ||. IX GRAVSAT - 1 GEOSTATIONARY (60° West of Center) 

A POSTERIOR! DENSITY SIGMAS 
21 SHORT DATA ARCS (DATA NOISE = ,05 CM/SEC) 

ADJUSTED PARAMETERS AND A PRIORI SIGMAS: UNADJUSTED PARAMETERS AND UNCERTAINTIES 

81 DENSITY BLOCKS, a = k MGAL lM» DENSITY BLOCKS, a = 4 MGAL 

LOW SATELLITE STATE, a [HCL]= (2M,5M,10M) 

HIGH SATELLITE STATE, a [HCL]= (20M, 50M, 100M) 
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TABLE III.X GRAVSAT - 1 GEOSTATIONARY ( 6 0° West of Center) 

A POSTERIORI DENSITY SIGMAS 
21 SHORT DATA ARCS (DATA NOISE = .05 CM/SEC) 

ADJUSTED PARAMETERS AND A PRIORI SIGMAS: UNADJUSTED PARAMETERS AND UNCERTAINTIES 

8] DENSITY BLOCKS, 0=4 MGAL 144 DENSITY BLOCKS, a = 4 MGAL 

LOW SATELLITE STATE, o [HCL.]= ( 1 OM, 25M.50M) 

HIGH SATELLITE STATE, a [HCL]= (j00M,250M,500M) 
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TABLE I II. XI GRAVSAT - 1 GEOSTATIONARY (60° West of Center) 

A POSTERIORI DENSITY SIGMAS 
21 SHORT DATA ARCS (DATA NOISE = .05 CM/SEC) 

ADJUSTED PARAMETERS AND A PRIORI SIGMAS: UNADJUSTED PARAMETERS AND UNCERTAINTIES 

81 DENSITY BLOCKS, a = b MGAL DENSITY BLOCKS, a = k MGAL 

LOW SATELLITE STATE, a [HCL]= (20M, 50M, 100M) 

HIGH SATELLITE STATE, a [HCL]= (200M.500M, 1 OOOM) 
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TABLE I I I. XI I 


GRAVSAT - 1 GEOSTATIONARY (60°West of Center) 

A POSTERIORI DENSITY SIGMAS 
21 SHORT DATA ARCS (DATA NOISE = .05 CM/SEC) 


ADJUSTED PARAMETERS AND A PRIORI SIGMAS: 

81 DENSITY BLOCKS, a = b MGAL 

LOW SATELLITE STATE, a [HCL]= {kOH, 100M.200M) 

HIGH SATELLITE STATE, a [HCL]= (400M, 1 000M,2000M) 


UNADJUSTED PARAMETERS AND UNCERTAINTIES 
144 DENSITY BLOCKS, a = k MGAL 






























































































































SET IV: 


Solution Set IV investigates the recovery of 2 1/2° x 2 1/2° density 
blocks with the "high-low" GRAVSAl'/GEOSTATIONARY configuration described 
in Section 6.1. Epoch state a priori uncertainties comparable to those of 
Tables II.I-II.VII are displayed in Tables IV.I-IV.V with a priori density 
block uncertainties of 4 mgals and 8 mgals. 

Table IV. I presents the most favorable case with virtually no error 
in the satellite epoch states and a 4 mgal a priori for the density block. 
Comparing with Table II. I shows a much poorer recovery ability for the 
2 1/2° x 2 1/2° resolution of the gravity fine structure. This result is 
more clearly demonstrated in Tables IV.II-IV.VI. Note that an assumed a 
priori density block uncertainty of 8 mgal and very accurate epoch state 
vectors (Table IV. V) results in a recovery error in the center adjusted 
block of 2.36 mgal, while halving the density block uncertainty to 4 mgal 
(Table IV. II) results in reducing the error to 1.77 mgal. Tables IV. VII 
and IV. VIII reproduce the error analysis study of Table IV. Ill (data 
noise .05 cm/sec) with the data noise assumed to be .01 cm/sec and .005 
cm/sec, respectively. This reduces the recovered uncertainty of 1.86 mgal 
to .89 mgal and .74 mgal, respectively. Table IV. IX displays the correlations 
of the adjusted blocks with the center block for the recovery depicted in 
Table IV. I. Relative to the correlations observed for 5° x 5° blocks of 
Sets I and II, the nearest neighbors in the north-south direction display 
a marked increase. However, these correlations indicate no difficulty in 
resolving blocks of 2 1/2° x 2 1/2° with the GRAVSAT at 250 km. altitude 
providing highly accurate SST data is available. Satellite-to-satellite tracking 
for the "high-low" GRAVSAT/GE0STATI0NARY configuration (with the GRAVSAT at 
250 km. altitude) does not significantly improve the gravitational fine 
structure at 2 1/2° x 2 1/2° resolution with data noise at .05 cm/sec. 

Until state-of-the-art instrumentation is available which will produce SST 
range-rate data at the .005 cm/sec noise level recovery of such resolution 
is not feasible. 
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TABLE IV. I GRAVSAT - GEOSTATIONARY (2.1/2° x 2 1/2° Blocks) 

A POSTERIORI DENSITY SIGMAS 
10 SHORT DATA ARCS (DATA NOISE = .05 CM/SEC) 

ADJUSTED PAREMETERS AND A PRIORI SIGMAS: UNADJUSTED PARAMETERS AND UNCERTAINTIES 

81 DENSITY BLOCKS, .a = A MGAL 1 DENSITY BLOCKS, a = A MGAL 

LOW SATELLITE STATE, a [HCL] = (.001M,.0plM, .001M) 

HIGH SATELLITE STATE, a [HCL] = (. 00 1M, . 001 M, . 001M) 







■ 

■ 

■ 

■ 

■ 

■ 

■ 

■ 

■ 

■ 






■ 

■ 

■ 

■ 

■ 

-.003 

■ 

■ 

■ 

■ 






■ 

■ 

■ 

■ 

■ 

-.012 

■ 

■ 

■ 

■ 






■ 


3.97 

3.88 

3.87 

3.87 

3.87 

ng 

3.96 

3.99 





.18 

.47 

.58 

.62 

.64 

m 

.27 

.12 



■ 

3.99 . 

3.90 

3.59 

.31 

3.58 

3.58 

3.58 

3.60 

3.88 

3.98 




.04 

.12 

.33 

.34 



.36 

.35 

.19 

.09 



■ 

3.98 

3.62 

fggj 


0ES 

Ppl 

T.S31 

3.64 

.19 

3.97 




.04 

.12 

■El 

.09 

.10 

. 2 2 1 

.08 



■ 

3.97 

BEEI 


mm 

1. 78 

1. 79 

2.12^ 

3.46 

BEEfi 




.05 

a 

dl 

a 

.05 

.04 

.171 

1 

.22 

.08 

-.003 

-.005 


3.96 

3.21 

.18 

*2.05 

m 

i . 

1.79 


wm 

RSI 

3.28 

3.95 

-.012 

-.007 

-.005 

.06 

.04 

■a 

. .04 

.15' 

| 

.23' 

.08 



■ 

3.96 

3.06 

mm 

1.93 

1.92 

1.96 

2. 23l 

3.15 

3.94 



■ 

.07 

.24 

m 

.04 

.06 

.05 

wm 

.28 

.09 



■ 

3.96 

3.14 

.25' 

m 

L1-. 

2.34 

2.32 

.ii 

2.36 

2.6lJ 

.-Hi 

3.19 

.27 

3.94 



■ 

.07 

.08 

.11 

.10 



■ 

3.98 

3.80 

3.63 

3.59 

.33 

3.59 

3.59 

3.64 

wm 

BBi 



■ 

.06 

.14 

.25 

.34 

.35 

.39 

■ 

mil 



■ 

3.99 

3.98 

3.92 

.34 

3.87 

mg 

3.87 

3.87 

3.92 

3.98 



■ 

.05 

.11 

.56 

.62 

.65 

.65 

.43 

.17 



■ 

■ 

■ 

■ 

■ 

m 

■ 

■ 

■ 

■ 



■ 



■ 

■ 

■ 

■ 

■ 

m 

■ 

■ 

■ 

■ 



■ 



■ 

■ 

■ 

■ 

■ 

.0006 

■ 

■ 

■ 

■ 



■ 


TABLE IV. II GRAVSAT - GEOSTATIONARY (2 1/2° x 2 1/2° Blocks) 

A POSTERIORI DENSITY SIGMAS 
10 SHORT DATA ARCS (DATA NOISE =• .05 CM/SEC) 

ADJUSTED PARAMETERS AND A PRIORI SIGMAS: UNADJUSTED PARAMETERS AND UNCERTAINTIES 

81 DENSITY BLOCKS, a = 4 MGAL 144 DENSITY BLOCKS, a = 4 MGAL 

LOW SATELLITE STATE, a [HCL]= (1M,2.5M,5M) 

HIGH SATELLITE STATE, a [HCL]= ( 1 OM, 25M, 50M) 
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TABLE IV. Ill GRAVSAT - GEOSTATIONARY (2 1/2° x 2 1/2° Blocks) 

A POSTERIORI DENSITY SIGMAS 
10 SHORT DATA ARCS (DATA NOISE = .05 CM/SEC) 

ADJUSTED PARAMETERS AND A PRIORI SIGMAS: UNADJUSTED PARAMETERS AND UNCERTAINTIES 

81 DENSITY BLOCKS, a = k MGAL ]kk DENSITY BLOCKS, c = b MGAL 

LOW SATELLITE STATE, a [HCL]® (2M,5M,10M) 

HIGH SATELLITE STATE, a [HCL]= (20M.50M, 100M) 
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TABLE IV. IV GRAVSAT - GEOSTATIONARY (2 1/2° x 2 1/2° Blocks) 

A POSTERIORI DENSITY SIGMAS 
10 SHORT DATA ARCS (DATA NOISE = .05 CM/SEC) 

ADJUSTED PARAMETERS AND A PRIORI SIGMAS: UNADJUSTED PARAMETERS AND UNCERTAINTIES 

81 DENSITY BLOCKS, a = A MGAL l^A DENSITY BLOCKS, a = A MGAL 

LOW SATELLITE STATE, a [HCL]= (AOM, 100M.200M) 

HIGH SATELLITE STATE, a [HCL]= (AOOM, 1 OOOM, 2000M) 
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TABLE IV. V GRAVSAT - GEOSTATIONARY (2 1/2° x 2 1/2° Blocks) 

A POSTERIORI DENSITY SIGMAS 
10 SHORT DATA ARCS (DATA NOISE = .05 CM/SEC) 

ADJUSTED PARAMETERS AND A PRIORI SIGMAS: UNADJUSTED PARAMETERS AND UNCERTAINTIES 

81 DENSITY BLOCKS, a = 8 MGAL 14A DENSITY BLOCKS, o = 8 MGAL 

LOW SATELLITE STATE, a [HCL]= (1M,2.5M,5M) 

HIGH SATELLITE STATE, a [HCL]= (10M,25M,50M) 

















































































































■ 

■ 

■ 

■ 

. 004 

B 

B 

B 

■ 


B 





■ 

■ 

■ 

■ 

.008 

■ 

a 

B 

■ 


1 





■ 

■ 

■ 

■ 

.016 

■ 

■ 

B 

■ 

! 






8.00 

7.98 

mm 

7.94 

7.94 

7.94 

7.95 

SB 

7.99 




.02 

.04 

.07 

.08 


.09 

.08 . 

II 

.03 




7.99 

.02 

7.85 

H9 

H 

7.00 


7.01 

BBS! 

7.83 

7.99 




.09 

.26 

.25 

.26 


.11 

.03 




7.99 

7.68 

5.73 

5 * 40 ~ 

574 0 ~ 
.09 

5742 

5 .*79 1 

7.68 

7.98 




.02 

.07 

.17 

■“ 

.10 

.23 

i 

.08 

.03 




wm 

6.94 

4.15 

3.79 

3.78 

3.81 

B 

7.01 

7.95 




E 

.13 

.15 

.09 

.10 

.10 

.17 

.05 

-.004 

-.006 

-.012 

7.95 

6.20 

4.09 

3 . 73 

a 

3.75 

4.18 1 

6.33 

7.93 

-% 018 ' 

-.010 


.04 

.20 

.10 

.06 

.07 

« 

.26 

m 





6.42 

4. 39 

3.98 

3.94 

4.00 


6.58 

7.95 




Id 

.14 

.15 

.09 

.10 

.11 

.20 1 
1 

.18 

.05 




7.97 

6.88 

EBiM 

5.77 

Hie 

5.81 

6.09 1 

6.91 

7.95 




.04 

.19 

m 

.10 

.11 

. 1.9 

.25 

.06 




wm 

7.58 

7.16 

7.04 

mm 

wm 

7.21 

nra 

7.96 




a 

.16 

.24 

.25 

.25 

.25 

.24 

.17 

.06 




7.99 

HU 

7.95 

7.94 

m 

7.94 

7.95 

HR! 

7.99 




.02 


.08 

.08 

BSE 

.09 

.08 

Bds 

.03 




■ 

■ 

■ 

■ 

.017 

■ 

■ 

■ 

■ 







■ 

■ 

■ 

■ 

.008 

■ 

■ 

■ 

■ 







■ 

■ 

■ 

■ 

.004 

■ 

■ 

■ 

■ 

1 




TABLE IV. VI GRAVSAT - GEOSTATIONARY (2 1/2° x 2 1/2° Blocks) 

A POSTERIORI DENSITY SIGMAS 
10 SHORT DATA ARCS (DATA NOISE = .05 CM/SEC) 

ADJUSTED PARAMETERS AND A PRIORI SIGMAS: UNADJUSTED PARAMETERS AND UNCERTAINTIES 

81 DENSITY BLOCKS, a = 8 MGAL 1 44 DENSITY BLOCKS, a = 8 MGAL 

LOW SATELLITE STATE, o [HCL]= (*»0M, 100M.200M) 

HIGH SATELLITE STATE, a [HCL]= (AOOM, 1 OOOM, 2000M) 
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TABLE IV. VI I GRAVSAT - GEOSTATIONARY (2 1/2° x 2 1/2° Blocks) 

A POSTERIORI DENSITY SIGMAS 
10 SHORT DATA ARCS (DATA NOISE = .01 CM/SEC) 

ADJUSTED PARAMETERS AND A PRIORI SIGMAS: UNADJUSTED PARAMETERS AND UNCERTAINTIES 

81 DENSITY BLOCKS, a = 4 MGAL 144 DENSITY BLOCKS, a = 4 MGAL 

LOW SATELLITE STATE, a [HCL]= (2M,5M,10M) 

HIGH SATELLITE STATE, o [HCL]= (20M,50M, 100M) 
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TABLE IV.VIII GRAVSAT - GEOSTATIONARY (2 1/2° x 2 l/2°- Blocks) 

A POSTERIORI DENSITY SIGMAS 
10 SHORT DATA ARCS (DATA NOISE = .005 CM/SEC) 

ADJUSTED PARAMETERS AND A PRIORI SIGMAS: UNADJUSTED PARAMETERS AND UNCERTAINTIES: 

81 DENSITY BLOCKS, a = b MGAL 1A4 DENSITY BLOCKS, a = b MGAL 

LOW SATELLITE STATE, a [HCL]= (2M,5M,10M) 

HIGH SATELLITE STATE, a [HCL]= (20M, 50M, 100M) 
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TABLE IV. IX GRAVSAT - GEOSTATIONARY (2 1/2° x 2 1/2° Blocks) 

•CORRELATIONS OF AOJUSTED BLOCKS WITH CENTER BLOCK 
10 SHORT DATA ARCS (DATA NOISE = .05 CM/SEC) 

ADJUSTED PARAMETERS AND A PRIORI SIGMAS: UNADJUSTED PARAMETERS AND UNCERTAINTIES: 

81 DENSITY BLOCKS, a «* k MGAL 1M DENSITY BLOCKS a = k MGAL 

LOW SATELLITE STATE, a [HCL]= (.001M, .001M, .001M) 

HIGH SATELLITE STATE, a [HCL]= (.001M, .001M, .001M) 




























































SET V: 


Solution Set V investigates the use of the "high-high-low" GRAVSAT/ 
GEOSTATIONARY configuration described in Section 6.1 for the recovery of 
high latitude 5 6 x 5° equal-area blocks. Tables V.I-V.VI repeat the a priori 
uncertainties of the epoch states and density blocks found in the studies of 
Tables II.I-II.VII and Tables III . I-III . VII . Making the appropriate corrections 
to the resultant a posteriori adjusted block errors by dividing by to account 
for the difference in data passes (see Section 6.1), we find gravity fine 
structure recovery comparable to the 5° x 5° equal-angular blocks located 
near the equator. 


SET VI: 


Solution Set VI as described in Section 6.1 is an extension of the "low- low" 
results obtained in report [A] and is included here for purposes of comparison 
with the "high-low" configurations of Sets I and II. Tables VI.I-VI.VII are 
designed to be compared with Tables I.I-I.VII and II.I-II.VII. The "high-low" 
results are superior in every case, although the degree of superiority is 
dependent on the assumed orbital knowledge of the low satellite. For the most 
favorable case assuming virtually error-free epoch states, the "high-low" mission 
is superior by a factor of approximately 20, while for assumed GRAVSAT epoch 
state uncertainties of (lm, 2.5m, 5m) and (2m, 5m, 10m) in HCL coordinates, 
the factor is approximately 10 and 7, respectively. For assumed GRAVSAT 
epoch state vectors with uncertainties of (40m, 100m, 200m) the superiority 
factor is only approximately 1.5. However, this comparison is somewhat 
misleadingly optimistic with regard to the "low-low" concept due to the fact 
that the orbital accuracy assumed is much more difficult to obtain for 
"low-low" short data arcs (see report [A]) than for "high-low" data arcs. 

Table VI. VIII displays the correlations for the error analysis study shown 
in Table VI. I, and reveals much higher values in the north-south direction 
(especially for nearest and next nearest neighbors) than the corresponding 
values for the "high- low" configurations of Sets I and II. 
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TABLE V.l GRAVSAT - 2 GEOSTATIONARY (POLAR 5* EQUAL-AREA BLOCKS) 

A POSTERIORI DENSITY SIGMAS 
12 SHORT DATA ARCS (DATA NOISE = .05 CM/SEC) 

ADJUSTED PARAMETERS AND A PRIORI SIGMAS: UNADJUSTED PARAMETERS AND UNCERTAINTIES: 

77 DENSITY BLOCKS, o = k MGAL 108 DENSITY BLOCKS, a = 4 MGAL 

LOW SATELLITE STATE, a [HCL]= ( . 001 M, . 001 M, . 001 M) 

HIGH SATELLITE STATE, a [HCL]= (.001M, .OOIM, . OOIM) 
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TABLE V. II GRAVSAT - 2 GEOSTATIONARY (POLAR 5° EQUAL-AREA BLOCKS) 

A POSTERIORI DENSITY SIGMAS 
12 SHORT DATA ARCS (DATA NOISE = .05 CM/SEC) 

ADJUSTED PARAMETERS AND A PRIORI SIGMAS: UNADJUSTED PARAMETERS AND UNCERTAINTIES 

77 DENSITY BLOCKS, ° = k MGAL 108 DENSITY BLOCKS, a = U MGAL 

LOW SATELLITE STATE, ° [HCL]= (1M,2.5M,5M) 

HIGH SATELLITE STATE, ° [HCL]= (1 0M,25M,50M) 
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TABLE V. I I I 


GRAVSAT - 2 GEOSTATIONARY (POLAR 5° EQUAL-AREA BLOCKS) 

A POSTERIORI DENSITY SIGMAS 
12 SHORT DATA ARCS (DATA NOISE = .05 CM/SEC) 


ADJUSTED PARAMETERS AND A PRIORI SIGMAS: 

77 DENSITY BLOCKS, a = b MGAL 

LOW SATELLITE STATE, o [HCL]= (2M,5M,10M) 

HIGH SATELLITE STATE, a [HCLJ= (20M.50M, 1 OOM) 


UNADJUSTED PARAMETERS AND UNCERTAINTIES: 
108 DENSITY BLOCKS, 0=4 MGAL 
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TABLE V.IV GRAVSAT - 2 GEOSTATIONARY (POLAR 5° EQUAL-AREA BLOCKS) 

A POSTERIORI DENSITY SIGMAS 
12 SHORT DATA ARCS (DATA NOISE = .05 CM/SEC) 

ADJUSTED PARAMETERS AND A PRIORI SIGMAS: UNADJUSTED PARAMETERS AND UNCERTAINTIES: 

77 DENSITY BLOCKS, <; = 4 MGAL 108 DENSITY BLOCKS, a = 4 MGAL 

LOW SATELLITE STATE, a faCL]= (10M.25M.5OM) 

HIGH SATELLITE STATE, 0 [HCL]= (1 00M.250M.500M) 
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TABLE V.V . GRAVSAT - 2 GEOSTATIONARY (POLAR 5° EQUAL-AREA BLOCKS) 

A POSTERIORI DENSITY SIGMAS 
12 SHORT DATA ARCS (DATA NOISE = .05 CM/SEC) 

ADJUSTED PARAMETERS AND A PRIORI SIGMAS: UNADJUSTED PARAMETERS AND UNCERTAINTIES: 

77 DENSITY BLOCKS, o = 4 MGAL 108 DENSITY BLOCKS, a = k MGAL 

LOW SATELLITE STATE, a [HCL]= (20M, 50M, 1 OOM) 

HIGH SATELLITE STATE, a [HCL]= (200M.500M, 1000M) 
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TABLE V.VI GRAVSAT - 2 GEOSTATIONARY (POLAR 5° EQUAL-AREA BLOCKS) 

A POSTERIORI DENSITY SIGMAS 
12 SHORT DATA ARCS (DATA NOISE = .05 CM/SEC) 

ADJUSTED PARAMETERS AND A PRIORI SIGMAS: UNADJUSTED PARAMETERS AND UNCERTAINTIES: 

77 DENSITY BLOCKS, a = 4 MGAL 1 08 DENSITY BLOCKS, a = k MGAL 

LOW SATELLITE STATE, o [HCL]= (40M, 1 00M,200M) 

HIGH SATELLITE STATE, o [HCL]= (400M, 1000M,2000M) 
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TABLE VI. I GRAVSAT - GRAVSAT (LOW-LOW) 

A POSTERIORI DENSITY SIGMAS 
21 SHORT DATA ARCS (DATA NOISE, = .05 CM/SEC) 

AND A PRIORI SIGMAS: UNADJUSTED PARAMETERS AND UNCERTAINTIES 

a = k MGAL 144 DENSITY BLOCKS, o = 4 MGAL 

a [HCL]= ( . 001 M, , 001 M , . 001 M) 


ADJUSTED PARAMETERS 
8l DENSITY BLOCKS, 
SATELLITE STATES, 
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TABLE VI. I I GRAVSAT - GRAVSAT (LOW-LOW) 

A POSTERIORI DENSITY SIGMAS 
21 SHORT DATA ARCS (DATA NOISE = .05 CM/SEC) 

ADJUSTED PARAMETERS AND A PRIORI SIGMAS: UNADJUSTED PARAMETERS AND UNCERTAINTIES 

81 DENSITY BLOCKS, a = 50 MGAL ]kk DENSITY BLOCKS, a = 50 MGAL 

SATELLITE STATES, a [HCL]= ( . 001 M, .001M, . 001M) 
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TABLE Vl.lll GRAVSAT - GRAVSAT (LOW-LOW) 

A POSTERIORI DENSITY SIGMAS 
21 SHORT DATA ARCS (DATA NOISE = .05 CM/SEC) 

ADJUSTED PARAMETERS AND A PRIORI SIGMAS: UNADJUSTED PARAMETERS AND UNCERTAINTIES 

81 DENSITY BLOCKS, a = k MGAL lV* DENSITY BLOCKS, o = k MGAL 

SATELLITE STATES, ' a [HCL]= (lM,2.5M,5M) 
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7.0 SUMMARY 


This study applied covariance error analysis techniques to investigate 
estimation strategies for the "high-low" mission for accurate local recovery 
of gravitational fine structure from multiple short data arcs, considering the 
effects of aliasing. The low drag-free GRAVSAT spacecraft was configured in 
a circular polar orbit at 250 km. altitude, while various configurations of 
the GEOPAUSE and geosynchronous spacecraft were utilized as high relay satellites. 
Sets of 5° x 5° equal-angular blocks near the equator and 5° x 5° equal-area 
blocks near the poles were used to represent the surface density block model 
of the high order geopotential used in the study. 

Satellite-to-satellite relative range-rate data was used in the error 
analysis investigation, as range data was found to be much less sensitive to 
fine structure recovery. A data noise of .05 cm/sec and a data rate of one 
observation every 5 seconds were assumed for most of the error analysis 
studies. No ground tracking data was included and satellite state accuracy 
requirements of the estimation of the density blocks were investigated in 
the error analysis by a priori error sigmas for the epoch states. The sen- 
sitivity of the data type to small a priori errors in the low satellite epoch 
state velocities led to severe aliasing errors in the estimation of local sets 
of density blocks, forcing the GRAVSAT epoch state parameters to be included 
with the surface density parameters in the estimation process. The estimation 
of local sets of density blocks was much less sensitive to the epoch state 
errors of the high relay satellite, so that for reasonable orbital uncertainties 
(order of tens of meters) the high satellite epoch state parameters could be 
treated as unadjusted in the estimation process. 

The error analysis for the recovery of local sets of density blocks using 
independent short data arcs demonstrated that the estimation strategy of 
simultaneously estimating a local set of density blocks covered by data and 
two surrounding "buffer layers" of density blocks not covered by data resulted 
in almost negligible aliasing errors on the estimation accuracy for blocks 
near the center of the recovered set due to unadjusted surface dnesity blocks. 
Estimates of blocks away from the center of the recovered set should be 
discarded. Error analysis studies using equal-area blocks at high latitudes 
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(near the pole) gave results comparable to those obtained with density blocks 
near the equator. An extrapolation to 5° x 5° resolution from a numerical 
simulation performed for global density block recovery using long data arcs 
indicates that local recoveries, using estimation strategies to reduce 
aliasing errors, can produce comparable recovery accuracies. 

Tables 7. 1-7. 3 present a comparison summary between the recovery accuracies 
for the "high-low" GRAVSAT/ GEOPAUSE and GRAVSAT/ GEOSTATIONARY configurations 
and the "low-low" configuration. The results presented in mgals are relative 
to gravity anomaly blocks. They are obtained by multiplying the results of 
Section 6.0 by 2ir (see Section 2.0). The summary tables clearly display 
the superiority of the "high-low" mission over the "low-low". For the 
low satellite assumed accurate to 2 meters radial, 5 meters cross-track, 
and 10 meters along-track, and the geostationary satellite assumed accurate 
to 20 meters radial, 50 meters cross-track and 100 meters along-track, then a 
priori density block errors of 25 mgals can be reduced to 2 mgals. For these 
orbital errors multiplied by 20, the a priori density block errors can be reduced 
to 8 mgals. 

These recovery accuracies are comparable to those obtained by Argentiero [1] 
in simulating gravity anomaly recovery using gradiometer data. We conclude, 
contrary to the findings for the "low-low" SST mission, that knowledge of the 
fine structure of the earth’s gravitational structure can be significantly 
improved at the 5° x 5° resolution using satellite-to-satellite tracking data 
from a "high-low" satellite configuration. However, reasonably accurate orbits 
must be available for the low GRAVSAT. 

The study demonstrates that the "high-low" SST data with the low. satellite 
at 250 km. altitude and a measurement noise of .05 cm/sec does not significantly 
improve the gravitational fine structure at the 2 1/2° x 2 1/2° resolution. A 
factor of ten improvement in the measurement noise (to .005 cm/sec) is needed 
to establish the feasibility of gravity fine structure recovery at 2 1/2° 
resolution from SST range-rate data. 
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TABLE 7.1 


GRAVSAT-GEOSTATIONARY 



STATE a* 

s 


GRAVITY a's 


(meters) 



(mgals) A; Priori = 25 mgals 


H 

C 

L 


LOW: 

1 

2.5 

5 

1.2 

HIGH: 

10 

25 

50 


LOW: 

2 

5 

10 

1.7 

HIGH: 

20 

50 

100 


LOW: 

10 

25 

50 

5.2 

HIGH: 

100 

250 

500 


LOW: 

20 

50 

100 

7.2 

HIGH: 

200 

500 

1000 


LOW: 

40 

100 

200 

8.5 

HIGH: 

400 

1000 

2000 
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TABLE 7.2 


GRAVSAT-GEOPAUSE 


STATE a's 

(meters) 

H C L 

LOW: 2 5 10 

HIGH: .412 


GRAVITY a's 

(mgals) A Priori = 25 meals 


1.88 


LOW: 

10 

25 

50 

HIGH: 

2 

5 

10 

LOW: 

20 

50 

100 

HIGH: 

4 

10 

20 

LOW: 

20 

50 

100 

HIGH: 

200 

500 

1000 

LOW: 

40 

100 

200 

HIGH : 

400 

1000 

2000 


5.4 


7.16 


7.48 


8.67 
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TABLE 7.3 


GRAVSAT- ("LOW-LOW") 



STATE o ’ 
(meters) 

s 


GRAVITY a ' s 

(mgals) A Priori = 25 mgals 


H 

C 

L 


LOW: 

1 

2.5 

5 

11 

LOW: 

2 

5 

10 

11.7 

LOW: 

10 

25 

50 

12.4 

LOW: 

20 

50 

100 

12.5 

LOW: 

40 

100 

200 

12.5 
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8.0 NEW TECHNOLOGY 


The effort under this contract consisted of the application of covariance 
error analysis techniques for investigation of estimation strategies for the 
"high-low" SST mission for accurate local recovery of gravitational fine 
structure, considering the effects of aliasing errors. Frequent reviews and 
a final survey for new technology were performed. It is believed that the 
mathematical and programming techniques and algorithms developed do not 
represent "reportable items", or patentable items, within the meaning of the 
New Technology Clause. Our reveiws and final survey found no other items 
which could be considered reportable items under the New Technology Clause. 
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APPENDIX A. GLOBAL LONG-ARC SOLUTIONS FOR THE "LOW-LOW" CONFIGURATION 


The results of report [A] demonstrated that aliasing from unadjusted 
density blocks led to severe corruption of the estimates of adjusted blocks 
when long data arcs were used in a local recovery. To establish the relation- 
ship between accuracies obtainable from global recoveries (no aliasing from 
unadjusted density blocks) and a global solution obtained as a composite of 
short-arc local recoveries, an error analysis experiment was performed to 
simulate a global 5° x 5° equal-area block recovery via long-arcs. 

The experiment was to scale the "low-low" satellite geometry and the 
rotation rate of the earth so as to obtain global error analysis solutions 
for smaller and smaller density block sizes and to observe the character of 
the global solutions as they approached the 5° x 5° block local solutions of 
report [A]. The satellite configuration consisted of two spacecraft in polar 
circular orbits separated by the length of the equatorial equal-area block 
and at altitudes equal to the length of the equatorial block. The rotation 
rate of the earth was set to 


2tt 
~ nP 

where n is the number of equal-area blocks on the equator and P is the 

satellite orbital period. The data rate was such that there were four measurements 

per pass per block at the equator, and a single data arc of duration 

T a = nP 
A 

was used which produced two passes at the equator per block. This experiment 
was performed for global sets of 90° (n=4) , 45° (n=8) , 30° (n=12) , 20° (n=18) , 

15° (n=24) and 10° (n=36) equal-area blocks. The SST relative range-rate 
measurements were given .05 cm/sec standard deviation, the a priori sigmas of 
the density blocks were set at 1 mgal and 50 mgal, and the satellite epoch 
states were included as adjusted parameters in the simulated recovery. 

Figures A. 1.1 - A. 1.3 display the a posteriori error sigmas for the most 
favorable case of 20° x 20° equal-area blocks for 1 mgal and 50 mgal a priori 
on the density blocks, respectively, with virtually perfectly known orbital 
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epoch states. The smaller errors near the poles reflect the more dense data 
coverage. Figure A. 1.3 displays the correlation pattern of the global 
set of 105 20° x 20° blocks with block number 44. 

Figure A. 1.4 displays the extrapolation of the sigmas for the equatorial 
blocks of all the solutions assuming perfectly known orbital states to a 
5° x 5° global set for 1 mgal and 50 mgal a priori. The extrapolated value 
of .8 mgal for the case of 1 mgal a priori compares favorably with the value 
for the center 5° x 5° equal-angular block of Table II. 3 of report [A] (.48 mgal) 
when corrected for the difference in data coverage between the two solutions, 
i. e. , 


Factor = 


' fp asses per block report [A] Data points per block report [A] 

\l Passes per block Appendix A Data points per block Appendix A 



» 2.74, 

so that the global recovery error becomes 

.8 mgal/FACTOR - .30 mgal. 

This result indicates the equivalence of local short-arc recoveries for gravi- 
tational fine structure using the "low-low" configuration, when strategies are 
employed to reduce aliasing, with global recoveries using long data arcs. 
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Adjusted Parameters and A Priori SigmasJ 

105 20° Equal -Area Blocks, o = 50 MGAL 
Satellite States, a[HCL] = (.001M, .001M,.001M). 






















A POSTERIORI SIGMAS (MGAl) 


I 


• A Priori Surface Density Slock Sigmas = 1 MGAL 
0A Priori Surface Density Block Sigmas = 50 MGAL 



, t 1 1 1 * 1 1 1 1 | i — 

0 .10 20 30 40 50 60 70 80 90 

EQUAL-AREA BLOCK SIZE 

FIGURE A. 1.4 Recovery Accuracy For Long Arc "Low-Low" Configuration With Global Coverage Over 
Equal-Area Density Blocks. 


